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Abstract 

The harmonic analysis of anisotropy of Ultra High Energy Cosmic Rays is 
performed for the Dark Matter halo model. In this model the relic superheavy 
particles comprise part of the Dark Matter and are concentrated in the Galac- 
tic halo. The Ultra High Energy Cosmic Rays are produced by the decays of 
these particles. Anisotropy is caused by the non-central position of the Sun 
in the Galactic halo. The calculated anisotropy is in reasonable agreement 
with the AGASA data. For more precise test of the model a comparison of 
fluxes in the directions of the Galactic Center and Anticenter is needed. 
The spectrum of Ultra High Energy Cosmic Rays (UHECR) is measured now up to a 
maximum energy of (2 — 3)10 20 eV P],0]. More than 1000 particles are detected at energies 
higher than 1 ■ 10 19 eV ||£|,fj,^]. The detailed energy spectrum was recently presented in 
||. No steepening of the spectrum has been observed in the energy range between 10 18 
and 2 • 10 20 eV. If extragalactic, the UHE protons must have the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff at energy E 1 / 2 = 6 • 10 19 eV 0. A similar cutoff should exist if primaries 
are extragalactic nuclei [j9|,j3[|10] or photons JXTJJT^ • 



Recently it was suggested that UHECR can be generated by the decay of superheavy 
relic particles [|T3|-|T6|. These particles can be effectively produced in the post-inflationary 



Universe ]T3],|n],[nj and can constitute now a small or large part of Cold Dark Matter (DM) . 
As any other form of Cold DM these relic particles are concentrated in the halo of our galaxy, 



and thus UHECR produced by their decays do not exhibit the GZK cutoff [13]. Realistic 
particle candidates for this scenario and possible mechanisms to provide the long lifetime 
for superheavy particles are discussed in Ref. Hl3 j |l4| ,|l9|,p0| . 



The halo model discussed above has three signatures 0,0,^31 : the excess of high energy 
photons in the primary radiation, direct signal from a nearby clump of DM (e.g. Virgo 
Cluster), and anisotropy caused by the asymmetric position of the Sun in the Galactic halo. 
These signatures allow to confirm or to reject the DM halo hypothesis by the data of existing 
arrays. 

As calculations show [p^ , [23|| , anisotropy reveals itself very strongly in the direction of 
the Galactic Center. This prediction can be reliably examined by the Pierre Auger detector 



in the southern hemisphere [23]. At present there is no detector which can observe this 
direction. In this article we present the calculations of anisotropy for the arrays in the 
northern hemisphere, taking as an example the geographical position of the Yakutsk and 
AGASA arrays. 

We shall assume that primary photons dominate in the decays of SH particles as QCD 
calculations [25] imply. Then the flux of UHE photons in the direction (£, 0) per unit solid 



angle can be written as 



/(C) = K / drp x (R), (1) 



n 



where r and R are the distances to a decaying X-particle from the Sun and the Galactic 
Center, respectively, ( is the angle between the line of observation and the direction to the 
Galactic Center, is the azimuthal angle in respect to Galactic plane (the flux depends 
only on (), px{R) is the mass density of superheavy particles (X) at a distance R from 



the Galactic Center, K is an overall constant, r max (() = J Rf t — r% sin £ + r cosC, r Q 



8.5 kpc is a distance between the Sun and Galactic Center, Rh is the size of the halo, in our 
calculations we shall use Rh = 50 kpc (the values of 100 and even 500 kpc result in similar 
anisotropy); the distance R from the Galactic Center to the decaying particle is given by 
R 2 (C) = r 2 + Tq — 2rr Q cos £. 

We shall use two distributions of DM in the halo: one given by the Isothermal Model 

(iso) a, 



and the other - following the NFW numerical simulation |27| 



P{R) ~ R/R s (l + R/R s r (3) 

In the ISO model we shall use for R c the values 5 kpc, 10 kpc and 50 kpc. For the NFW 
model the calculations are performed for R s equal to 30 kpc, 45 kpc and 100 kpc. The NFW 



distribution |27| is given in terms of the virial radius r 20 o, the rotational velocity at the 
virial distance V200, the constant 5 C and the dimensionless Hubble constant h. We applied 
this distribution to our Galaxy using the following parameters: the local density of DM 
Pdm(j q ) = 0.3 GeV/cm 3 , V200 = 200 km/ s and h = 0.6. As a result we obtain R s w 45 kpc. 

The flux ([]]) was expressed first in terms of galactic coordinates, longitude / = and 
latitude b, which is given by cos b = cos£/ cos0, and then transferred into equatorial co- 
ordinates, declination 5 and right ascension a. We calculated the amplitudes of the first 
and the second harmonics (Ai and A 2 , respectively) and the phase a of the first harmonic 
for the geographical position of the Yakutsk and AGASA arrays. These quantities are the 
standard ones used for measured anisotropy. The results are given in Table 1 (predictions 
for the AGASA array are shown in brackets). Depending on the parameters of the DM dis- 
tribution, the anisotropy varies from 10% to 45%. The phase of the first harmonic a ~ 250° 
is close to the RA of Galactic Center, a ~ 265°. 

After this paper was submitted for publication we receved the preprint by Medina- Tanco 



and Watson , where similar calculations were performed. The results of both calculations 



are displayed in Fig.l for E > 4-10 19 eV together with the data of AGASA (AG) and Yakutsk 
( YK) arrays. The AGASA anisotropy is taken from analysis of Ref . [E3] . Our calculations 



(BM) agree well with that of Ref. [ffijfl . Both agree with the data of AGASA array and do 
not contradict to the Yakutsk data. 
Acknowledgements 

The authors are grateful to P.Blasi, M.Hillas, B.Hnatyk, M.Nagano, P.Sokolsky, A.Vilenkin, 
and A.Watson for discussions and correspondence. 



4 



REFERENCES 

N.Hayashida et al, Phys. Rev. Lett. 73 (1994) 3491. 
[2] D.J.Bird et al., Ap.J., 424 (1994) 491. 
[3] G.Cunningham et al, Ap.J. 236 (1980) L 71. 
[4] S.Yoshida et al., Astroparticle Physics 3 (1995) 105. 

[5] B.N.Afanasiev et al., in Extremely High Energy Cosmic Rays: Astrophysics and Future 
Observations , edited by M.Nagano (ICRR, University of Tokyo, 1996) p. 32. 

[6] M.Takeda et al, Phys.Rev.Lett. 81 (1998) 1163. 

[7] K.Greisen, Phys. Rev. Lett. 16 (1966) 748; 

G.T.Zatsepin and V.A.Kuzmin, Pisma Zh. Exp. Teor. Fiz. 4 (1996) 114 . 

V.S.Berezinsky, S.V.Bulanov, V.A. Dogiel, V.L.Ginzburg, and V.S.Ptuskin, "Astro- 
physics of Cosmic Rays", chapter 4, ELSEVIER, 1990. 

[9] V.Berezinsky, S.Grigorieva and G.T.Zatsepin, Proc. 14th Int. Cosm. Ray Conf. (Munich 
1975), v.2, p.711. 



[10] L.N.Epele and E. Roulet, J.High Energy Phys. 9810 (1998) 009, |astro-ph/9808i~04 

[11] V.S.Berezinsky, Soviet Phys. Nucl. Phys. 11 (1970) 399 . 

[12] R.J.Protheroe and P.L.Biermann, Astroparticle Physics 6 (1996) 45. 

[13] V.A.Kuzmin and V.A.Rubakov , Talk at the Workshop "Beyond the Desert", Castle 
Rindberg 1997; Yadernaya Fisika 61, 1122 (1998). 

[14] V.Berezinsky, M.Kachelriess and A.Vilenkin, Phys. Rev. Lett. 79 (1997) 4302. 

[15] P.H.Frampton, B.Keszthelyi and N.J.Ng, |astro-ph/ 9709080. 



[16] M.Birkel and S.Sarkar, Astrop. Phys. 9 (1998) 297. 



[17] D.J. H.Chung, E.W.Kolb and A.Riotto, Phys.Rev. D 59 (1999) 023501. 

[18] V.Kuzmin and I.Tkachev, JETP Lett. 68 (1998) 271. 

[19] K.Benakli, J.Ellis and D.Nanopoulos, hep/ph/9803333. 

[20] K.Hamaguchi, Y. Nomura and T.Yanagida, Phys.Rev. D 58 (1998) 103503 

[21] S.L.Dubovsky and P.G.Tinyakov, JETP Lett. 68 (1998) 107. 

[22] V.Berezinsky, PBlasi and A.Vilenkin, Phys. Rev. 58 (1998) 103515. 

[23] A.Benson, A.Smialkowski and A.W.Wolfendale 1998, submitted to Astrop. Phys. 

[24] J.W.Cronin, Nucl. Phys. B (Proc. Suppl.) 28B, 213, (1992). 

[25] V.Berezinsky and M.Kachelriess, Phys. Lett. B 434 (1998) 61. 



[26] A.V.Kravtsov, A.K.Klypin, J.S.Bullock and J.R.Primack. |astro-ph/9708176| (to be pub- 
lished in Ap.J). 

[27] J.F.Navarro, C.S.Frenk and S. D.M.White, Astroph. J. 462, 563 (1996); 
F.F.Navarro |astro-ph/980TU75. 



[28] G. A. Medina- Tanco and A.A.Watson, 1998, submitted to Astrop. Phys. 



6 



TABLES 



TABLE I. Anisotropy 



ISO model NFW model 



R c 


A x 


A 2 


a 


R s 


Ax 


A 2 


a 


5 kpc 


0.43(0.46) 


0.11(0.13) 


250° 


30 kpc 


0.41(0.45) 


0.10 (0.13) 


250° 


10 kpc 


0.32(0.35) 


0.06(0.07) 


250° 


45 kpc 


0.37(0.41) 


0.09 (0.11) 


250° 


50 kpc 


0.15(0.15) 


0.01(0.01) 


250° 


100 kpc 


0.33(0.36) 


0.07 (0.08) 


250° 



FIGURE CAPTIONS 

Figure 1: Amplitude and phase of the first harmonic of anisotropy for the AG AS A 
and Yakutsk arrays. Solid lines are for the ISO distribution of DM and dots (NFW) - for 
the NFW numerical simulation. BM and TW show the results of this paper and Ref. ]28 



respectively. The three dots of BM (from left to right) are given for R s = 30, 45 and 
100 kpc, respectively; five dots of TW - for R s = 10, 20, 30, 50 and 100 kpc. The 



AGASA data are taken from Ref. PS . 
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